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ABSTRACT 

We use a sample of 42 supernovae detected with the Advanced Camera for Surveys on-board the 
Hubble Space Telescope as part of the Great Observatories Origins Deep Survey to measure the rate of 
core collapse supernovae to z ~ 0.7 and type la supernovae to z ~ 1.6. This significantly increases the 
redshift range where supernova rates have been estimated from observations. 

The rate of core collapse supernovae can be used as an independent probe of the cosmic star formation 
rate. Based on the observations of 17 core collapse supernovae, we measure an increase in the core 
collapse supernova rate by a factor of ~ 1.6 in the range 0.3 < z < 0.7, and an overall increase by a 
factor of ~ 7 to z ~ 0.7 in comparison to the local core collapse supernova rate. The increase in the rate 
in this redshift range in consistent with recent measurements of the star formation rate derived from 
UV-luminosity densities and IR datasets. 

Based on 25 type la supernovae, we find a SN la rate that is a factor 3 — 5 higher at z ~ 1 compared 
to earlier estimates at lower redshifts (z < 0.5), implying that the type la supernova rate traces a higher 
star formation rate at redshifts z > 1 compared to low redshift. At higher redshift (z ^ 1), we find 
a suggested decrease in the type la rate with redshift. This evolution of the la rate with redshift is 
consistent with a type la progenitor model where there is a substantial delay between the formation of 
the progenitor star and the explosion of the supernova. Assuming that the type la progenitor stars have 
initial main sequence masses 3M < M < 8M , we find that 5 — 7% of the available progenitors explode 
as type la supernovae. 



Subject headings: 



galaxies: distances and redshifts — galaxies: stellar content 
- surveys — 



supernovae: general 



1. INTRODUCTION 

The rates of supernovae (SNe) at different redshifts pro- 
vide important information about the evolution of a num- 
ber of physical processes over cosmic times. As core col- 
lapse SNe (i.e., type II and Ibc supernovae, hereafter CC 
SNe) originate from massive short-lived stars, the rate of 
these events should reflect on-going star formation and 
therefore offers an independent way to determine the cos- 
mic star formation rate (SFR). Furthermore, the SN rate 
(SNR) directly probes the metal production at different 
cosmological epochs. 

Type la SNe do not directly follow the SFR since there 
is a delay between the formation of the progenitor star and 
the explosion of the SN. Constraining this delay time is im- 
portant for a better understanding of the processes leading 
to the type la explosions, and therefore essential for the 
usefulness of type la SNe as cosmological distance indica- 
tors. With accurate measurements of the SN la rate, it 
should be possible to set constraints on this poorly known 
delay time as discussed in e.g., Madau, Delia Valle, & 
Panagia (1998) and Dahlen & Fransson (1999). A first 



attempt to constrain the delay time using observations is 
presented in Gal- Yam & Maoz (2004) and Maoz & Gal- 
Yam (2004). 

Using SNRs as evolutionary probes requires the rates to 
be measured over cosmological distances where the prop- 
erties of the Universe, e.g., the SFR(z), are believed to 
change by a significant fraction. In practice, this means 
to redshifts z ~1, which is higher than has previously 
been feasible. Local and nearby SN la rates have been 
measured in the field by numerous surveys, e.g., at red- 
shifts < z >~ 0.01 (Cappellaro, Evans, & Turatto 1999), 

< z >= 0.1 (Hardin et al. 2000), < z >= 0.11 (Strolger 
2003), and < z >= 0.114 (Reiss 2000). Similarly, at larger 
distances, rates have been measured at < z >= 0.38 (Pain 
et al. 1996), < z >= 0.46 (Tonry et al. 2003), and 

< z >= 0.55 (Pain et al. 2002). While rates at z £ 0.4 are 
in general higher than the local rates, these measurements 
are still consistent with a constant type la rate out to 
z ~0.6 (Tonry et al. 2003). In clusters of galaxies, rates 
have been estimated at < z >— 0.25 and < z >= 0.9 (Gal- 
Yam, Maoz, & Sharon 2002). 
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By contrast, CC SNe are typically ~ 2 mag fainter and 
searches are affected by severe selection biases, thus only 
locally determined rates of CC SNe exist (e.g., Cappellaro 
et al. 1999; Strolger 2003). 

The Great Observatories Origins Deep Survey (GOODS, 
Giavalisco et al. 2004a) has offered an unprecedented op- 
portunity to obtain a sample of high rcdshift SNe, ideal for 
studying the evolution of the SNR with rcdshift. Using the 
Advanced Camera for Surveys (ACS) on-board the Hub- 
ble Space Telescope (HST), GOODS has detected 42 SNe 
in the rcdshift range 0.2 ;$ z ;$ 1.6. Repeated observations 
of two fields, combined with the high spatial resolution of 
the ACS, have resulted in a sample with well understood 
detection efficiency and systematics. In this paper, we de- 
rive and discuss the rates of CC and type la SNe based on 
this unique sample. In an accompanying paper (Strolger et 
al. 2004), we describe details of the search, including data 
reduction procedures, SN finding and classification, selec- 
tion biases, detection efficiency and completeness. This 
paper also presents a detailed investigation of SN la pro- 
genitor delay time distributions derived from the GOODS 
SN sample. The SN sample, including magnitudes and 
redshifts, is presented in Strolger et al. and Riess et al. 
(2004). Constraints on cosmological parameters using the 
GOODS SN sample are presented in Ricss et al. 

The data are presented in § 2. In § 3 we describe how 
the SN rates are determined from observations. Results 
are given in § 4, followed by an investigation of the rela- 
tion between star formation and SN rates in § 5. Results 
are discussed in § 6. Finally, conclusions and a summary 
are given in § 7. 

Throughout this paper, we adopt a flat lambda dom- 
inated cosmology (ft\ — 0.7, 17m = 0.3) and a Hubble 
constant Hq = 70 km s _1 Mpc -1 . Magnitudes are given 
in the Vega based system. 

2. DATA 
2.1. Observations 

During the GOODS campaign, two fields (HDF-N and 
CDF-S) were observed using HST/ACS during five epochs 
separated by ~45 days. Both fields were observed with 
multiple pointings, covering an effective area ~150 sq. ar- 
cmin per field. The search was conducted using observa- 
tions in the z-band (F850LP filter), with additional obser- 
vations in V and i (F606W, F775W), allowing colors to 
be calculated. For a number of high rcdshift SN la candi- 
dates, near IR photometric (HST/NICMOS) and spectro- 
scopic (VLT, Keck) follow up was conducted, as reported 
in Riess et al. (2004) and Strolger et al. (2004) 

When searching for SNe, we subtracted two consecu- 
tive epochs creating a residual image, which thereafter was 
both automatically and visually scanned for objects. Each 
area was independently investigated by two teams. Com- 
pleteness was estimated by adding artificial point sources 
with different magnitudes to the real images. The result- 
ing completeness limit is rriF850LP — 25.9, defined as the 
magnitude where 50% of these sources where recovered. 

2.2. The Supernova Sample 

The GOODS SN sample consists of 42 SNe detected in 
both CDF-S and HDF-N. For 29 of these, we have spectro- 



scopic redshift determinations from either the SNe them- 
selves, or the host galaxies, while we have photometric 
redshifts (Mobasher et al. 2004) for an additional 12 SNe. 
We only lack redshift for one SN with no apparent host 
galaxy. 

We wish to differentiate between type la and CC SNe 
when calculating rates since these types originate from dif- 
ferent progenitors and explosion mechanisms. It is there- 
fore essential to conduct accurate type determination for 
the full sample. We divide the SNe into types using all 
available information including spectra, redshifts of SNe or 
host galaxies, magnitudes, colors and light curve shapes. 
We find that 25 SNe are consistent with being type la, 
with the remaining 17 being CC SNe. We are most con- 
fident in the type determination of 20 of the la and 7 of 
the CC SNe. For the remaining ones, the type determina- 
tion is less certain and a few may be misclassified. When 
calculating rates, we account for this by including possible 
misclassifications in the systematic errors. 

3. DETERMINING SN RATES FROM OBSERVATIONS 

In this paper, we express SNRs in units of number of 
exploding SNe per rest-frame year and comoving volume 
element (i.e., yr _1 Mpc -3 ). Alternatively, rates can be ex- 
pressed in supernova units, SNu's (1 SNu = 1 SN per 100 
year per 10 10 Lq), which are normalized to the i?-band 
luminosity density, and are preferentially used for local 
measurements. We do not use SNu's because this intro- 
duces further uncertainty in surveys covering cosmological 
distances where the rest-frame B-band luminosity and its 
evolution has to be determined to provide proper rates 
in SNu's. When we compare local rates given in SNu's 
with rates per volume, we convert units assuming a local 
B-band luminosity density 2.0xl0 8 ftX Q Mpc~ 3 , approxi- 
mately evolving as (1 + z) 19 at z <1 (Strolger 2003). The 
rcdshift correction is, however, small since the highest red- 
shift we convert is z ~ 0.1. 

We use Monte Carlo (MC) simulations to calculate the 
underlying SNR that is consistent with observations. This 
is done separately for CC and type la SNe. We start by as- 
suming an input SNR and use this to calculate the number 
of SNe exploding each year within a volume given by the 
field-of-view and the redshift range of interest. Each SN is 
assigned a peak magnitude, a random epoch on the light 
curve, a host galaxy extinction and a redshift. The ap- 
parent magnitude is thereafter calculated at this epoch on 
the light curve, as well as at an earlier epoch separated in 
rest-frame by 45/(1 + z) days, corresponding to the 45 day 
spacing between observations. By subtracting the fluxes 
between the two epochs, we derive a detection magnitude, 
Am, and the direction of the magnitude change, i.e., if 
the SN is rising or declining. This gives us a list of simu- 
lated SNe with known detection magnitudes and redshifts. 
The MC simulation is repeated 10,000 times to calculate 
the mean number of detectable SNe for a specified search 
set-up. Finally, to determine the observed rate, we adjust 
the input rate to match the number of detected SNe in a 
chosen rcdshift bin. 

Calculating the apparent magnitude of the SNe in the 
simulations includes specifying a number of SN character- 
istics. After the following general formula, we specify the 
CC and SN la characteristics in the two next subsections. 
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The apparent magnitude in the observed filter F850LP 
for a SN at rcdshift z which is at time t relative to the 
light curve peak in the observers frame is given by 

m F850L p(t, z) = M peaKX + AM x [t x (1 + z)" 1 ] 

+ DM(z) - K F850LP [z,tX (1 + z)- 1 ] + A A (1) 

where M peak ^\ is the peak absolute magnitude in a rest- 
frame filter A, AM\ is the light curve decline relative to the 
peak, DM(z) is the distance modulus. The K-correction, 
Kp S5aLP , is calculated using the formalism in Kim, Goo- 
bar, & Pcrlmuttcr (1996) and includes information on the 
spectral SED of the SN together with appropriate filter 
transmission curves and QE for the detector used. Fi- 
nally, A\ is the extinction in the SN host galaxy. This is 
further described below. 

3.1. Core Collapse SN Characteristics 

CC SNe can be divided into a number of different sub- 
types, characterized by e.g., light curve decline, spectral 
features or peak magnitude. We adopt the division into 
subtypes given in Richardson et al. (2002). Here CC SNe 
are divided into Ibc, IIL, IIP and Iln, where the Ibc and 
IIL types are further divided into a normal and a bright 
population, while the IIP type is divided into a normal 
and a faint population. A summary of peak i?-band mag- 
nitudes and corresponding dispersions are given in Table 1. 
We use the intrinsic fractions of subtypes given in Dahlen 
& Fransson (1999), with the addition of 1% of the SNe as- 
signed to a bright population as proposed by Richardson 
et al. This bright population is further divided so that half 
are Ibc and half IIL. Note that the division into different 
types is highly uncertain. Also, the intrinsic fractions may 
change with redshift (metallicity) . In § 5.1.1, we investi- 
gate how the subdivision into subtypes affects the results. 

To obtain K-corrections, we use the approach in Dahlen 
& Fransson (1999), where the SN spectra are charac- 
terized by black bodies with temperatures evolving from 
~25,000K near explosion to ~5,000K at late stages. We 
use these models since there does not exist a full sample of 
observed CC SN spectra of different subtypes and epochs 
on the light curve covering a sufficiently large range in 
wavelength to calculate K-corrections. For particular SNe 
(e.g., 1999em) with good observational coverage, we find 
that the blackbody approximation and observed spectra 
yield similar K-corrections, with deviations mostly < 0.1 
mag. This is typically less than the photometric uncer- 
tainty in the observed CC SNe, justifying the use of the 
blackbody approximation. Light curves for Ibc, IIL, and 
IIP are taken from Filippenko (1997), while for Iln, we use 
a light curve intermediate between IIL and IIP. 

CC SNe should preferentially occur in star bursting re- 
gions within the host galaxy. We therefore use a starburst 
extinction law (Calzetti et al. 2000) when calculating ab- 
sorption of the SN due to dust in the host galaxy. Further- 
more, we use a mean color excess E(B — V)=0.15. The de- 
tails of dust extinction, especially in star forming regions, 
are, however, highly uncertain and in § 5.1.1, we thor- 
oughly investigate how different assumptions about the 
extinction affects results. 

Each SN is placed in a disk with random inclination. 
The extinction is thereafter calculated by integrating the 



light path through the host galaxy using the recipe in 
Hatano, Branch, & Deaton (1998). This leads to a distri- 
bution of extinctions that peaks at low values, Ab < 0.1, 
but that also shows a tail with high extinction values 
(Ab > 5 — 10), originating from highly inclined host 
galaxies. 

3.2. Type la Characteristics 

To characterize type la SNe, we use the peak magnitude 
distributions in Tonry et al. (2003). Here type la SNe are 
divided into three Gaussians with peak magnitudes Mb=- 
19.6, -19.3, -17.8 and dispersions ct=0.30, 0.45, and 0.50, 
respectively (see Table 1). The SNe are distributed within 
these three distributions so that 20% are 'bright', 64% are 
'normal', and the remaining 16% are 'dim' (Li et al. 2001). 

K-corrections are calculated using a set of observed spec- 
tral templates from SN1994D, observed at days -5, 0, 5, 9, 
17, and 55 relative to peak in light curve. We interpolate 
between these spectra to get K-corrections at intermediate 
epochs on the light curve. We use a mean i?-band light 
curve taken from Riess et al. (1999). 

Absorption in the host galaxy is calculated as in Hatano 
et al. (1998). We divide the SN la so that 25% arc placed 
in a bulge, while the remaining are placed in a disk as 
suggested by Tonry et al. (2003). This is consistent with 
the fraction ^30% found by Farrah et al. (2002), as well 
as preliminary results from the GOODS SN sample. The 
scale height used for the disk in Hatano et al. (1998) is 
larger than for the CC SNe, leading to lower mean ab- 
sorption for SN la. For both bulge and disk components, 
the distribution of host galaxy B-band absorption has a 
very dominating peak at low extinctions (Ab < 0.1). 
There is, however, also a narrow tail out to high extinc- 
tions (Ab Si 4), especially for SNe occurring in inclined 
disks. The overall distribution of extinctions can be ap- 
proximated by an exponential, P(Ab) oc e~ AB , consistent 
with the observed distribution of SN la absorption in Jha 
et al. (1999). 

4. RESULTS 

4.1. Observed Core Collapse Rates 

Using the method described above, we calculate the CC 
SNR in two rcdshift bins 0.1 < z < 0.5 and 0.5 < z < 0.9, 
which is the first measurement of CC SN rates at cosmolog- 
ical distances. The resulting rates are 2.51 *q 75 -1 86 anc ^ 
3.96 tlm -2.6O in thc two bins > respectively (Table 2). 
Rates are given in units of 10~ 4 yr -1 Mpc~ 3 ft,i! . First er- 
rors quoted represent 68.3% confidence intervals, and are 
calculated from a distribution of rates based on 10,000 
MC simulations. The 95% confidence interval is well rep- 
resented by two times the quoted statistical errors. The 
second quoted errors are systematic and include e.g., the 
possibility of misclassification of SN type. Systematic er- 
rors are further discussed in §5.1.1. 

The measured rates are higher by factors ~ 4 and ~ 7 
at < z >= 0.3 and < z >= 0.7 compared to the local 
(z ~ 0.0) rate ~ 0.59 in Cappellaro et al. (1999) (we 
convert the local rate from supernova units (SNu's) as de- 
scribed in § 3). The rate in Cappellaro et al. is calculated 
using an empirical correction for inclination effects biasing 
against finding SNe with high extinction. Cappellaro et al. 
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compare this method with the corrections in Hatano et al. 
(i.e., similar to what we use) and find results that are con- 
sistent within ~15%. The main result, that our rates at 
0.3 < z < 0.7 are significantly higher than the local rate, 
therefore seems robust, as further discussed in §5.1. 

With the CC SNe being directly related to the SFR, the 
increase in the SNR reflects an increase in star formation 
over this redshift range, as discussed below. 

4.2. Observed Type la rates 

We calculate the type la SNR in four redshift bins 
0.2 < z < 0.6, 0.6 < z < 1.0, 1.0 < z < 1.4, and 
1.4 < z < 1.8, using the MC simulations described 
above. Results are given in Table 2. We find a rate 
0.69 toM io.25 at < z >= 0.40 (rates given in units 
10 _4 yr _1 Mpc _3 /i7 ). This rate is somewhat higher than 
rates previously measured at z ~ 0.4-0.5 (Pain et al. 1996, 
2002; Tonry et al. 2003). The differences are within the 
error-bars, however, it may be possible that the derived 
ground-based rates are underestimated due to systematic 
effects, as further discussed in § 5.2. At < z >— 0.80, 
we derive a rate 1.57 j^' 2 5 ~-o 53' significantly higher than 
previous measurements at lower redshift. At even higher 
redshifts, we find a rate that is consistent with declin- 
ing at z > 1. The rate at < z >= 1.20 is, however, still 
higher than local measurements. In numbers, we find rates 



1.15 



+0.47 +0.32 
-0.26 -0.44 



and 0.44 ±° f 5 t%\i at < z >= 1.20 and 
< z >= 1.60, respectively. Errors represent statistical and 
systematic errors, as described in §4.1. 

While earlier measurements were still consistent with a 
type la SNR constant with redshift, our new results sug- 
gest that the SNR increases by a factor ~ 5 from the local 
Universe to a cosmic time corresponding to z ~1, i.e., 
when the Universe was less than half its present age, and 
thereafter the rate shows a slight decrease toward higher 
redshift. Fitting a constant type la rate to the combined 
data results in a reduced chi-squarc x 2 /^=2.9, equivalent 
to rejecting a constant rate with ~99.9% probability. The 
derived chi-square is based on statistical errors only, tak- 
ing systematic errors into account should lower the signif- 
icance. 

The increased rate to redshift z ~1 is expected, since it 
reflects the observed higher SFR at redshifts z £ 1 — 2, 
compared to local rates. Nevertheless, this is the first 
time it has been measured, supporting the general shape 
of the star formation history as predicted from e.g., mea- 
surements of galaxy UV-luminosities. We revisit this dis- 
cussion below. 

5. RELATION BETWEEN STAR FORMATION AND 
SUPERNOVA RATES 

With sufficient knowledge of how a star becomes a SN, 
it is possible to calculate the relation between the SFR 
and the SNR. Depending on the nature of the relation, it 
is possible to either use the SNR as a probe for the SFR, 
or to use independently determined SFRs and SNRs to 
investigate the relation between the two and extract infor- 
mation about the nature of the SN progenitors. For CC 
SNe, the relation between SFR and SNR is fairly straight- 
forward since the progenitors are massive stars with life 
times spanning only a small fraction of the Hubble time. 
Therefore, if the mass range of CC SN progenitors and the 



IMF arc known, we can calculate a constant that relates 
the SFR to the SNR. For type la SNe, the case is more 
complicated since there is a more substantial delay time 
involved, representing the time elapsed between the for- 
mation of the progenitor star and the explosion of the SN 
la. The distribution of delay times is not well understood, 
but is expected to be a large enough fraction of the Hubble 
time to untie any direct relation between SFR and SNR. 
Also, the efficiency, i.e., fraction of available progenitor 
white dwarfs that actually explode as SN la, is unknown. 
In this paper, we use the recently determined SFR in 
Giavalisco et al. (2004b) as our fiducial model (Ml) when 
investigating the relation between observed SNR and un- 
derlying SFR. In Giavalisco et al., the SFR at z > 3.5 is 
determined using GOODS/ACS data, while rates from the 
literature are used at lower redshift (Lilly et al. 1996; Con- 
nolly et al. 1997; Stcidcl et al. 1999). The SFR, which 
is derived from UV-luminosity densities, after correcting 
for dust extinction, is characterized by a factor ~ 10 in- 
crease between z = and z ~ 2, thereafter the rate de- 
clines slowly toward higher redshifts up to z ~ 6. The 
SFR fitted to these data is shown in Figure 8 in Strolger 
et al. (2004). For comparison, we also fit data points in 
Giavalisco et al. without corrections for dust extinction. 
This model (M2) shows a similar increase in the SFR to 
z ~ 1.5 as in model Ml. At higher redshifts, model M2 
shows a decrease in the SFR, similar to the shape of the 
SFR first suggested by Madau et al. (1996). 

5.1. Star Formation Rates from, Supernova Rates 

The progenitors of CC SNe arc believed to be massive 
stars in the range 8M < M < 50M Q (e.g., Nomoto 1984; 
Tsujimoto et al. 1997). For an assumed initial mass func- 
tion (IMF), -0(M), the number of CC progenitors per unit 
formed stellar mass is 
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/X HM)dM 
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Mi>(M)dM 



(2) 



Evaluating equation (2), assuming a Salpeter IMF 
(Salpeter 1955), yields fc=0.0069 M' 1 . Since the life time 
of the progenitor star is sufficiently short, we can assume 
a direct relation between the mass of formed stars (SFR in 
units M Q yr~ 1 Mpc -3 ), and the number of exploding CC 
SNe (SNR in units yr _1 Mpc -3 ), 



SNR(z) = kxh 2 x SFR(z) 



(3) 



In this paper, we compare our estimated SNRs with SFRs 
derived from observed luminosity densities. Calculating 
the SFR from luminosity densities results in a rate propor- 
tional to the Hubble constant to the power one, i.e. oc h. 
However, measuring SNRs from counting SNe in a volume, 
results in a rate that is proportional to h 3 . Therefore, rates 
measured in these two ways will differ by a factor /i 2 , as 
shown in equation (3). 

In Figure 1, we plot our measured rates at < z >= 
0.3 and < z >— 0.7 as filled circles, together with the lo- 
cal rate at z ~ 0.0 from Cappellaro et al. (1999) (filled 
square). The SNR in Figure 1 is given by the scale on 
the left y-axis, while the corresponding SFR is given by 
the right y-axis. The relative off-set between the scales is 
given by the factor k x h 2 in equation (3). Also plotted 
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are the two different SFR models described above. The 
SFR derived from the CC SNe is slightly higher, but still 
consistent with the SFR derived from extinction corrected 
UV-luminosities shown as model Ml. Based on our two 
data points, we find that the SFR increases by a factor 
~ 1.6 between < z >= 0.3 and < z >— 0.7. However, 
a flat rate over this range is consistent within error-bars. 
Including the local rate from Cappcllaro ct al. (1999), we 
find that the SN based SFR shows a steep (factor ~ 7) in- 
crease in the redshift range < z ^ 0.7, similar to model 
Ml. In more detail, a chi-square fit shows that the SFR 
derived from SNe is on average ~ 13% higher than the SFR 
in model Ml, this difference is, however, within errors and 
may not be significant. Shown in the figure is also the dust 
enshrouded SFR derived from analysis of various mid- and 
far- infrared datasets (Chary & Elbaz 2001). We hereafter 
refer to this as the IR SFR. The IR SFR matches the SN 
derived SFR at z ~ 0.3, but is a factor ~ 2 higher at 
z ~ 0.7. Considering errors, we conclude that both the IR 
and the extinction corrected UV-luminosity derived SFRs 
are consistent with our data points. We also note that if 
we combine the local rate from Cappellaro et al. (1999) 
with our data points, we can reject a flat CC SNR with 
99.9% probability (\ 2 jv — 7.3), considering statistical er- 
rors only. The conclusion that the rate increase with red- 
shift, therefore seems robust. We further investigate this 
conclusion in §5.1.1, where systematic errors discussed. 

For comparison, we also plot the SNR assuming no dust 
extinction (Figure 1, open circles, see also Table 2). This 
is certainly not physically correct since we know that the 
light from CC SNe is at least partly extinguished, but it 
shows an absolute lower limit to the rates. Despite the 
absence of correction, the SNR we measure is still higher 
than the corrected local rate. The SFR derived from the 
uncorrected SNR approximately follows the shape of the 
uncorrected SFR model M2, but is on average a factor 
~1.5 higher. Finally, we also plot in Figure 1 the SFR 
derived by Lanzetta et al. (2002). We have here taken 
their 'middle' model and transformed rates to the cosmol- 
ogy used here. No correction for extinction is made in 
Lanzetta et al., we therefore compare their rates with the 
SFR derived from the uncorrected SNR. We find that the 
rates in Lanzetta ct al. are a factor ~ 2 lower than the 
rates we find, this difference is at a ~ 3ct level (99.6%) 
considering statistical errors only. There arc alternative 
models in Lanzetta et al. with somewhat higher SFR, 
however, all models predict rates lower than found here. 

5.1.1. Systematic Effects in CC SN Rates 

When calculating the CC and type la SN rates our re- 
sults depend on a correct type determination of our SN 
sample. As mentioned in § 2.2, we are most confident in 
the type determination of 7 out of 17 of the CC SNe and 
20 out of 25 of the la SNe (called 'Gold' and 'Silver' SNe 
in Strolger et al. 2004). To investigate possible system- 
atic effects that misclassification may introduce, we calcu- 
late errors in the rates assuming that all SNe that do not 
have absolutely certain determination of type have been 
assigned a wrong type (i.e., 'Bronze' type SNe in Strolger 
et al.). The systematic errors due to possible misclassifi- 
cation are 2.51 ii'Jg and 3.96 iS'Ji in the two bins, re- 
spectively (units are 10~ 4 yr~ 1 Mpc _3 ft,7 ). Inspecting the 



errors shows that statistical errors most often dominate 
over these systematic errors. There are, however, cases 
where the latter dominate. Note that the systematic er- 
rors due to possible misclassification are not lcr-errors, but 
represent a worst case scenario, i.e., where all SNe that we 
are not certain of, are misclassificd. In reality, we expect 
that at most a few SNe are misclassificd. 

Furthermore, when estimating the SNR, as well as a 
SFR from the SNR, our calculations include a number of 
assumptions and SN properties that, to some extent, are 
poorly known and may introduce additional systematic 
errors. First, the SNR depends on a number of char- 
acteristics, including light curves, peak magnitudes, K- 
corrections and dust extinction. Second, when deriving 
the SFR from the SNR, we have to specify the constant k, 
as well as the Hubble constant, in equation (3). 

To determine how sensitive our results are to the as- 
sumptions that go into the calculations, we investigate 
how much we have to vary the assumptions in order for 
the results to change by more than la of the statistical 
errors (a s tat)- For CC SNe, ±la sta t corresponds to ±30% 
in the estimated rates (using approximately the mean of 
the statistical errors of the two data points). After this 
investigation, we make an estimate of the total systematic 
error. We first investigate possible sources for systematic 
errors in the derived CC SNR. 

• Peak magnitudes. CC SN peak magnitude distribu- 
tions are subject to some uncertainty because deriv- 
ing these from observations depend on assumptions 
about extinction (as well as observational uncer- 
tainties). The typical uncertainty in peak magni- 
tudes for different CC subtypes in Richardson et al. 
(2002) is 0.1-0.2 mag. In order to change the es- 
timated number of CC SNe with ±a s t a t, we have 
to change the assumed peak magnitudes of all sub- 
types, and in the same direction, by +0.3/-0.6 mag, 
which is unlikely considering the typical peak mag- 
nitude errors quoted. To further investigate this, we 
rerun our MC simulations incorporating the peak 
magnitude uncertainties for each subtype given in 
Richardson et al. We find that this only introduces 
an extra dispersion in the resulting rates by 2 and 
4% in the low and high redshift bins, respectively. 
This is insignificant compared to statistical errors. 

• Subtypes. Changing the division into different sub- 
types may affect results. It is, however, not simple 
to quantify this in a direct way. First, we exam- 
ine how the results change under the assumption 
that all CC SNe are of a single type represented 
by the SN IIP characteristics. We find that the re- 
sulting number of expected SN decreases by ~15%, 
which is less than the statistical errors and suggests 
that results are not highly dependent on the exact 
division into subtypes. The reason for this weak de- 
pendence is that omitting both the bright and the 
faint sub-populations cancels most of the net effect. 
What could change the results systematically, is if 
the fractions of very faint or very bright subtypes 
are significantly wrong. To examine the effects of 
overestimating either of these populations we cal- 
culate the predicted number of SNe after setting 
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the fractions of bright and faint subtypes to zero, 
respectively. Excluding the bright Ibc and IIL pop- 
ulation only changes the result by 3%, while an ex- 
clusion of the faint IIP population changes the rate 
by at most 15%, both changes are less than statis- 
tical errors. Next we examine how much we have 
to increase the fractions of faint or bright SNe in 
order to change the predicted number by a s t a t ■ We 
find that the intrinsic fraction of bright SNe have 
to increase from 1% to ~25%, while the faint frac- 
tion have to increase from 15% to ^35%. Neither 
of these scenarios seem likely. 

• Dust extinction. The resulting rates are dependent 
on the assumptions about dust extinction, as shown 
in Figure 1 where our extinction corrected rates are 
a factor ~ 2 higher than the uncorrected. In or- 
der to change our results by ±a s t a t, we have to 
change the assumed mean E(B — V)=0.15 with 
±AE(B — F)=0.06. Using a Galactic extinction 
law (Ry — 3.1) instead of a Starburst (Ry = 4.05) 
decreases the estimated rates by ~ 9-15%, signifi- 
cantly less than statistical errors. 

Besides these possible systematic effects when deriving the 
SFR, there are also systematic effects that may enter the 
derivation of the SFR from the SNR due to uncertainties 
in k and h in equation (3). 

• IMF. The constant k in equation (2) depends on the 
IMF and the assumed mass range for the progenitor 
stars. Changing the shape of the IMF assuming a 
Scalo IMF (Scalo 1986) instead of the Salpeter IMF, 
decreases k by a factor 2.6. However, changing the 
IMF also changes the conversion from observed UV- 
lurninosities to estimated SFRs. For a Scalo IMF, 
the conversion factor increases by a factor ~2, i.e., 
the net effect of changing the IMF is therefore a 
decrease in k by ^30%, which increases the mea- 
sured SFR by the same amount. Again this does 
not exceed the statistical error in the derived rates. 

• Progenitor mass range. We have chosen 50Mq as 
the upper limit to the progenitor mass since it is 
believed that more massive star become black holes 
without exploding as SNe (Tsujimoto et al. 1997). 
Due to the steep slope of the IMF, increasing the 
upper limit to 125M Q , only increases the factor k by 
7%. If we instead lower the upper limit to 30M Q , 
we get a resulting decrease in k by 9%. A more sig- 
nificant change in k comes from changing the lower 
mass limit of the progenitor stars. To change k cor- 
responding to a s t a t, we have to increase the lower 
integration limit to 10.2M Q . This could be possible 
since the lower mass range of CC SNe is not well 
constrained and is estimated to be in the interval 
8 - 11M Q (Timmes, Woosley, & Weaver 1996). 

• Hubble constant. Finally, the relation between SNR 
and SFR depends on the Hubble constant to the 
second power. According to Spergel et al. (2003), 
the error in the determination of Hg is less than 5%, 
suggesting that the uncertainty due to the h 2 fac- 
tor in equation (3), should be ~10%, clearly less 



than statistical errors. Note that the comparison 
between the SFRs derived from luminosity densi- 
ties and SNRs respectively, provides a method for 
determining Ho- For a stringent determination of 
Ho, both statistical and systematic errors must, 
however, be significantly better constrained than is 
possible today. At face values, we can take our ob- 
served rates and compare to the SFR derived from 
UV-luminosity densities to calculate Ho. This ex- 
ercise results in Ho=66±8 km s^ 1 Mpc -1 , where 
errors are statistical and represent la. 

Errors introduced by systematic effects described above 
are mostly non-Gaussian, and it is difficult to estimate the 
total effect of the uncertainty on the results. For the SNR, 
only assumptions about dust extinction may have a signif- 
icant effect on the rates. When deriving the SFR from 
the SNR, we add sources with possible systematic effects. 
To estimate the total error, we make the assumption that 
the uncertainty in peak magnitudes introduces a 5% er- 
ror, while division into subtypes introduce an additional 
15% uncertainty. Furthermore, we assume uncertainties 
due to dust extinction 10% and 20% in the two redshift 
bins, respectively. This results in a total error in the SNR 
of - 19 - 25%, equivalent to - 0.6 - 0.8cr stat . For the SFR 
derived from the CC SNR, we add three possible sources of 
statistical errors as shown by the listed points above. As- 
suming that each of these introduce a 10-15% uncertainty 
would lead to a systematic error that is of the same order 
as the statistical error (~ l<r s tat)- 

To investigate if the strong significance of the increased 
SNR in the redshift range < z < 0.7 found in §5.1 de- 
creases when taking systematic effects into account, we as- 
sume systematic errors ~ 0.6 — 0.8a sta t, as derived above. 
In section §5.1, we found that a flat CC rate could be re- 
jected with 99.9% probability, based on statistical errors 
only. Adding the estimated systematic errors decreases 
this probability to 99.4%, which is still significant. 

We also note that the systematic effects discussed above 
are mostly independent of redshift. This implies that the 
relative increase in SNR over the redshift range observed 
should remain the same, and therefore supporting an in- 
crease in SFR, even if the absolute normalization of our 
determined rates may have a systematic offset. Only the 
assumptions about the amount of dust extinction have a 
strong redshift dependent effect on the derived rates. As- 
suming a higher extinction results in a steeper increase 
with redshift of the estimated rates, as can be seen in Fig- 
ure 1 by comparing the corrected and uncorrected data 
points. 

In summary, there are a number of possible sources for 
systematic errors in our estimates, however, most errors 
are relatively small, not exceeding the statistical errors. 
For the SNR, we estimate that the summed systematic er- 
rors should be smaller than the statistical, while the sys- 
tematic errors may be of the same order as the statistical 
when it comes to the SFR derived from the SNR. Since we 
have shown that systematic errors are unlikely to dominate 
over the statistical, and that they are mostly independent 
of redshift, we are confident that the increase in supernova 
and star formation rates we observe are true features. The 
main concern is the amount of dust extinction. The effects 
of changes in the dust extinction are further discussed in 
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§6. When we derive the SFR from the SNR, we also note 
that the direction of the errors is mainly to increase the 
observed rates, e.g., this is the case if the lower limit for 
CC progenitor mass is larger than 8M Q , if the IMF is 
changed to become less top heavy or if the amount of dust 
extinction is underestimated. 

We finally note that the quoted systematic errors in 
Table 2 are sums of the 19-25% uncertainty derived 
above, and the uncertainty due to possible misclassifica- 
tion. These errors are therefore non-Gaussian. 

5.2. Type la supernova rates 

Even though the physics behind type la SNe has been 
extensively investigated using both observations and theo- 
retical simulations, there is still a lack of understanding of 
the mechanisms that proceed the explosion of this SN type 
(see Livio 2000, 2001 for reviews). The evolution of the 
type la SNR should follow the SFR, but shifted towards 
lower redshift after taking a delay time into account. For 
a particular distribution of delay times, 3>(£d), where td is 
the time elapsed between the formation of the progenitor 
star and the explosion of the SN la, the SNR is given by 
a convolution of the SFR over delay times 



SNR Ia {t) = v [ SFR(t r ) x $(t - t')dt' 

Jtw 



(4) 



where t is the age of the Universe and tF is the time 
corresponding to the redshift, zf, where the first stars 
formed. In this paper we set zf—W. is is the number of 
SNe per unit formed stellar mass (Mq 1 ). In Strolger et 
al. (2004), we calculate predicted redshift distribution as- 
suming two different functional forms for the distribution 
of delay times. The first model has an e-folding delay time 
distribution, $(td) oc exp(—td/r), while the second model 
has a Gaussian distribution, 3>(£d) <* exp{ — [(td — T )/2c] 2 }. 
We call t the characteristic delay time. Note that this pa- 
rameter has different meanings in the two functions. The 
e-folding distribution always has the highest probability 
for td=0, and has ~ 1/3 of the SNe with td > r. The Gaus- 
sian distribution has the highest probability for td — T, and 
1/2 of the SNe have td > t. We use two different Gaussian 
models, one 'narrow' with a = 0.2t and one 'wide' with 
a = 0.5t. By comparing the observed redshift distribution 
with the predicted, we investigate in Strolger et al. (2004) 
which shapes of the delay times function, as well as which 
ranges of r are consistent with the observed distribution. 
We find that all models favor characteristic time scales, 
r, around or greater than ~ 3 Gyr. Models with delay 
times t < 2 Gyr can be rejected with 95% confidence. 
These constraints are similar to the ones found by Gal- 
Yam & Maoz (2004), who estimate the delay time using 
an e-folding distribution only. Using a different approach, 
Maoz & Gal- Yam (2004) argue that the delay-time must 
be short if the iron mass observed in clusters of galaxies 
originate from SN la. Further comparisons between these 
investigations are presented in Strolger at el. (2004). 

In Figure 2, we plot observed data points from this inves- 
tigation together with rates from Cappcllaro et al. (1999), 
Hardin et al. (2000), Strolger (2003), Reiss 2000, Pain 
et al. (1996, 2002) and Tonry et al. (2003). We have 
converted the first four rates from SNu's using a _B-band 
luminosity density given in § 3. We also show SNRs de- 
rived from equation (4) using the three different delay time 



models. For the 'narrow' and 'wide' Gaussians, we use the 
best-fitting value r = 4.0 Gyr, derived from fitting the 
models to the GOODS dataset in Strolger et al. (2004). 
For the e-folding we choose r = 5.0 Gyr (there is no best- 
fitting value for this functional form, but values below ^2 
Gyr are rejected with 95% confidence). We further use 
SFR model Ml in equation (4). In Strolger at al., we find 
that the 'narrow' Gaussian delay time distribution fits the 
GOODS data better than the other two models. From Fig- 
ure 2, we see that this fit is also consistent with the low 
redshift (z ;$ 0.1) data taken from the literature. Only 
the two data points at z ~0.5 deviate significantly from 
the fit and are about a factor ~ 2 lower than the model 
suggests. A possible explanation for this deviation is that 
ground-based searches may suffer from unaccounted in- 
completeness. B. Barris (2003, private communication), 
through a rcanalysis of light curves of all variable objects, 
has found significant incompleteness in the ground based 
2001 survey for SN la reported by Barris et al. (2004). 
As many as half of the bona fide SN la at z < 0.6 were 
not recognized as such during the campaign. The salient 
characteristic of many of these missed SN la is extremely 
close proximity to a relatively bright galaxy which causes 
them to be identified as potential AGN, making them dif- 
ficult to observe spectroscopically regardless of their true 
nature. It is not possible to know how much of this bias 
was present in the preceding Fall 1999 campaign reported 
in Tonry et al. (2003), but it indicates that the rate es- 
timates contained therein could be low by up to nearly a 
factor of two, which would bring the rate at z = 0.5 into 
much better agreement with that found here. 

The reduced chi-squares when fitting the GOODS data 
points to models are x 2 /i^=0.46, 1.01 and 1.46, for the 
'narrow' Gaussian, 'wide' Gaussian and c-folding distribu- 
tions, respectively The 'narrow' Gaussian has the low- 
est value since it best reproduces both the sharp increase 
at z < 1 and the decline at z > 1. The reduced chi- 
squares for the same distributions, including also rates 
from literature are x 2 /i^=7.4, 7.9 and 10.6, respectively. 
The relatively large chi-squares are mostly due to the de- 
viant points at z <~ 0.5. Excluding these points results in 
X /i/=l.l, 3.2, and 6.3, respectively. A similar improve- 
ment is also achieved if the z ^ 0.5 points are multiplied by 
a factor ~ 2. These results again favor the 'narrow' Gaus- 
sian as the best model. In Strolger et al., we concluded 
that delay time models predicting a large fraction of SN 
la with delay times shorter than ^2 Gyr are inconsistent 
with data. Any such model does not follow the decline at 
z > 1 that our observations suggest. 



5.2.1. Systematic effects in la rates 

The measured type la rate is basically affected by the 
same systematic effects as the CC rate. First we have pos- 
sible errors due to misclassification of SN type. These 



errors alone are 0.69 + n'-M, 1-57 



+0.47 
-0.26' 



1.15 



+0.08 
-0.20 



and 



0.44 1q oi' m t ne f our redshift bins, respectively (in units 
10 _4 yr -1 Mpc~ 3 ft,7 ). For the remaining sources of errors, 
there are a few factors that make the SN la rates less sub- 
ject to systematics compared to CC SNe. Characteristics 
like peak magnitudes, light curve shape and spectra are 
more universal, as well as better known, for la SNe, and 
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should not introduce any significant systematics. The ex- 
tinction corrections could introduce systematics that sig- 
nificantly affect the derived rates. However, the depen- 
dence on assumed extinction is likely to be smaller than 
for CC SNe since the type la SNe on average occur in en- 
vironments that are less affected by dust extinction. To 
examine the sensitivity to systematic effects, we examine 
the type la SNe estimates in a similar way to the CC SNe. 
For the type la SNe we investigate how much we have 
to change the assumptions that go into our calculation in 
order to change the estimated rates by ±lo~ s t a t- 

• Peak magnitude. To change the estimated rates in 
each of the four redshift bins with +l(J s tat, the peak 
magnitudes have to be more than 0.5 mag fainter 
than the values given in Table 3. Alternatively, the 
peak magnitude have to be more than one magni- 
tude brighter to change the estimated number by 
— lcstot in each bin. These required shifts are sig- 
nificantly larger than the accuracy with which the 
(mean) peak magnitude of la SNe is known. In this 
investigation, we do not include the observed effect 
that brighter SNe la have a slower evolution on the 
light curve, while fainter SNe evolves faster. To 
a first approximation, these two effects cancel out, 
leaving no net effect on the derived rates. However, 
to investigate if there is a net effect at high redshifts 
where only the brightest SNe are expected to be de- 
tected, we run a set of MC simulations accounting 
for this. We use the relation between peak mag- 
nitude and light curve stretch given in Pcrlmuttcr 
et al. (1997). Repeating the simulations above, we 
find than the derived rates in the four redshift bins 
increase by less than 4%, with the smallest increase 
in the highest redshift bin (1%). Therefore, we find 
that the peak magnitude light curve stretch relation 
has little effect on the derived rates. 

• Subtypes. Setting the fraction of cither the faint or 
bright SN la subtypes to zero does not change the 
estimated rates by more than ~ 10%. To change 
the estimated rates by \<J s tat, we have to increase 
the fraction of faint SN la from 16% to - 50%. Set- 
ting the fraction to 100% for the bright population 
changes the rates by less than \a stat . These calcu- 
lations show that the division into subtypes should 
not introduce systematics that are comparable to 
statistical errors. 

• K- corrections. We use the spectrum of SN1994D 
to calculate K-corrections since this is the only ob- 
served SN with sufficient spectral coverage avail- 
able over a large number of epochs on the light 
curve. The good UV-coverage of SN1994D is es- 
sential for deriving the K-correction to redshifts 
z ~ 1.8. SN1994 is, however, about ~0.3 mag bluer 
in (U — B) color than the mean (U — B) color in Lei- 
bundgut et al. (1988). To investigate if the unusu- 
ally blue color of SN1994D affects results, we rerun 
our MC simulations with a magnitude correction to 
the K-corrections. We start to add the correction at 
the redshift where the effective wave-length of the 
observed filter (z-band) reaches the effective wave- 
length of redshifted rest-frame B-band (i.e., z ~ 1). 



The correction increases linearly to reach 0.3 mag 
at the redshift where the observed band probes the 
effective wave-length of the redshifted U-h&nd (i.e., 
z ~ 1.6). At higher redshift, the correction is set to 
a constant 0.3 mag. We find that the effect of this 
correction on results is small. The derived rates in- 
crease by less than 5% in the two highest redshift 
bins when including this color correction (bins at 
z < 1 arc not affected at all). We therefore con- 
clude that the use of the unusually blue SN1994D 
when deriving K-correction should not affect results 
more than marginally. 

• Dust extinction. Using only half the nominal 
amount of dust extinction reduces the estimated 
rates in the four redshift bins by 5, 10, 14 and 19%, 
respectively. All these changes are less than sta- 
tistical errors. To increase the estimated rates by 
lc s tat, we need to increase the dust extinction by 
a factor ~ 3, compared to the extinction model in 
Hatano et al. (1998), used here. 

Based on the discussion above, we find that the system- 
atic errors in the type la rates should not exceed statisti- 
cal errors. Making the simple assumption that the three 
first listed sources introduces uncertainties ^10%, ~5%, 
and ~5%, while the uncertainty due to dust extinction 
increases from 10% in the lowest redshift bin to 20% in 
the highest redshift bin, results in a total added error 
~ 16 — 23%, which is less than statistical errors in all bins, 
and is on average less that 0.6cr stat . We are therefore con- 
fident that general evolution in the la rate we observe is a 
true feature. Quoted systematic errors in Table 2 include 
these uncertainties as well as uncertainties due to possible 
misclassification of SN type. 

5.3. White dwarf explosion efficiency 

Fitting the predicted model distributions of SN la to the 
observed sample includes determining the normalization, 
v 1 in equation (4). This number tells us how many type 
la SNe explode per unit formed stellar mass. Previously 
in this investigation, we have made no assumptions about 
the mass range of type la progenitor stars. If we make 
an assumption constraining the progenitor mass range, we 
can calculate the fraction of stars in this mass range that 
subsequently explode as SN la, which we here call the ef- 
ficiency n, via the relationship 






(5) 



The mass range of progenitor stars is set to 3M Q < M < 
8M Q (Nomoto et al. 1994). By fitting each of the three 
delay time models, <&(td), in equation (4), we estimate 
i/=1.0xl0~ 3 , 1.2xl0~ 3 , and 1.3xl0~ 3 for the 'narrow' 
Gaussian, the 'wide' Gaussian and the e-folding distribu- 
tions, respectively. Evaluating equation (5) results in an 
efficiency of r\= 4.9%, 5.6%, and 6.3%, with which progen- 
itor white dwarfs explode as SN la, depending on the delay 
time model. Furthermore, if we assume that all stars with 
main sequence masses 0.8M Q < M < 8M Q become white 
dwarfs, then our observations suggests that 0.6 — 0.8% of 
the total number of white dwarfs will explode as SNe la. 
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These numbers are consistent (within observational errors 
and the uncertainties associated with the progenitor mod- 
els) with population synthesis calculations (e.g., Yungcl- 
son & Livio 2000, and Yungclson, private communication), 
which predict this fraction to be ~0.5%. 

5.4. Ratio of the Core Collapse to Type la SN Rates 
In Figure 3, we plot the redshift evolution of the intrin- 



sic ratio between CC and SN la rates, r, 



(CC/Ia)- 



For CC, 



we use rates calculated by fitting SFR model Ml to our 
observed data points. SN la rates based on the three de- 
lay time models, each having a characteristic delay time t 
quoted in § 5.2. For the 'narrow' and 'wide' Gaussian, we 
get local ratios between CC and la rates rrcc/ia) ~ 4.0 
and 2.9, respectively. This is consistent with the local 
ratio ^3.5 derived by Madau et al. (1998) using a compi- 
lation of locally determined rates. For the e-folding delay 
time models we get a local ratio 2.2. This lower value is 
mostly an effect of the higher type la rates at low redshift 
predicted in this model. The CC to la ratio stays fairly 
constant at a value ~ 3 to z ~ 0.7. The ratio increases 
rapidly at higher redshifts. The reason for this behav- 
ior is that the relatively long delay time shifts the SN la 
to lower redshifts compared to the SFR (and hence CC 
SNR). Therefore, SN la becomes relatively more common 
at lower redshift, lowering the CC to la ratio. In the case of 
a short delay time, the la SNR and the SFR should closely 
follow each other, giving a ratio ri cc / Ia \ that would be rel- 
atively constant with redshift. Using MC simulations, we 
find that the statistical uncertainty in the derived ratio 
is ~ ±30%. Systematic errors add to these, however, the 
systematic errors are mostly independent of redshift and 
should therefore only change the normalization of the ratio 
and not the evolution with redshift. Therefore, the result 
that the ratio is fairly constant to z ~ 0.7 and thereafter 
increases, should be robust. Note that these results re- 
fer to intrinsic ratios of exploding SNc. In a magnitude 
limited search, the detected ratio will be lower since CC 
SNe are typically fainter and are on average more severely 
affected by dust extinction. 

6. DISCUSSION 

We have derived a SFR to z ~ 0.7 from measurements 
of the CC SNR that is ~ 13% higher than (but consistent 
with) the SFR based on evolution of extinction corrected 
UV- luminosity densities. There are claims that a large 
fraction of the cosmic SFR is hidden from measurements 
based on rest-frame UV-luminosities. Driver (1999) ar- 
gues that extreme selection effects bias against finding low 
surface brightness galaxies at z > 0.2. These systems 
may therefore contain a large amount of undetected star 
formation. SNe occurring in low surface brightness envi- 
ronments are readily detected. In fact, SNe in these envi- 
ronments should be easier to detect since there is no noise 
added due to a background source when the host galaxy 
is undetected. However, we do not find more than one CC 
SN (out of 17) that lacks an apparent host galaxy. This 
suggests that only a limited amount of star formation may 
be hidden in such environments. The lack of host-less SNe 
also suggests that there is not a large amount of star for- 
mation occurring in faint low mass systems at high redshift 
as is the case locally (Cowie et al. 1996). 



It is difficult to estimate corrections due to dust extinc- 
tion in the host galaxies. Because SNe are detected in the 
optical, instead of rest-frame UV, they are less affected by 
dust extinction, however, corrections are still sufficiently 
large to change the estimated rates considerably, as shown 
in § 5. Observations of the luminosity density at longer 
wavelengths are less affected by dust extinction and may 
therefore offer a way to more directly determine the star 
formation. Chary & Elbaz (2001) derive the SFR from 
mid- and far-IR observations and find a rate that follows 
the overall shape of the UV-determined rate, with a sharp 
increase toz ~ 1.5, and thereafter a constant or declining 
rate, similar to our model M2. The rate derived from the 
IR is, however, in general significantly higher than UV de- 
rived rate, a difference that increases with redshift. The 
IR derived SFR is a factor ~ 2 higher at z = 0.5 and a 
factor ~ 3 higher at z = 1, as compared to the extinction 
corrected UV-derived SFR from Giavalisco et al. (2004b). 
Our results from CC SNc match the IR SFR at z = 0.3, 
but is a factor ^2 lower at z = 0.7. The reduced chi- 
squares when fitting the UV as well as IR derived SFRs 
to the observed points are 0.5 and 16, respectively, imply- 
ing that the UV-derived SFR fits the observations better. 
However, an increased amount of dust extinction in the 
host galaxy will make observations more consistent with 
the IR data. To investigate this, we adjust the dust extinc- 
tion so that the SFR derived from CC SNe fits the IR SFR. 
We find that a mean extinction E(B — U)=0.40 (instead 
of the nominal E(B — V)=0.15), results in the best fit be- 
tween our SN derived SFR and the IR SFR. The reduced 
chi-square using this extinction is 3.5. 

So far, we have assumed a constant E(B — V) over 
the redshift range probed. It is, however, possible that 
E(B — V) evolves so that extinction increases with red- 
shift. This was suggested by Totani & Kobayashi (1999) 
who derived an evolution in Ay of ~ 0.2 mag between 
z ~ and z ~ 1 when integrated over all galaxy types. 
Fitting our two data points independently results in best- 
fitting values E{B-V)=Q.l& at z=0.3 and £(B-U)=0.50 
at z=0.7, which corresponds to an evolution in Ay of ~ 1.4 
mag in this redshift range. This is higher than the evolu- 
tion found by Totani & Kobayashi, but is more consistent 
with the analysis of Chary (2004) who has re-calculated 
the average extinction as a function of redshift based on 
the most accurate estimates of the luminosity density at 
far-infrared and ultraviolet wavelengths and finds that the 
best-fit E(B - V) evolves from 0.19 at z=0.3 to 0.3 at 
2=0.7. 

Finally, a large fraction of star formation may be hid- 
den in extremely extinguished environments such as cores 
of starbursts or ULIRGs (Blain ct al. 1999). It is hard to 
use optically detected SNe to probe this activity, since the 
absorption in these systems is typically Ay ~ 10 (Man- 
nucci et al. 2003). Therefore, SFRs based on both UV- 
luminosities and optically detected SNc, will miss most of 
this population. Targeted observations, searching for SNe 
in the NIR, may be the way to detect a possible abun- 
dant population of CC SNe in these objects (Mattila & 
Meikle 2001). Note that for type la SNe it is possible that 
the delay-time makes detection possible if sufficient time 
is available for the dust to settle or disperse before the 
explosion of the SN. 
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We have found good agreement between the SFRs de- 
rived from SNe and galaxy luminosity densities. We do, 
however, note that there may be corrections to the derived 
rates. With the systematic effects being of the same order 
of magnitude as the statistical errors, it is necessary to in- 
crease the statistical sample, as well as learn more about 
systematic effects, in order to significantly reduce errors. 
Note, however, that uncertainties involved in deriving the 
SFR from high redshift galaxy UV-luminosities are also 
most certainly affected by systematic effects and unknown 
details regarding e.g., dust extinction and IMF, leading to 
uncertainties that should be of the same order of magni- 
tude. Therefore, measuring the SFR from observations of 
CC SNe, provides an important independent probe for the 
cosmic SFR. Furthermore, compared to measurements of 
the SFR based on UV-luminosities, the point source SNe 
are not affected by surface brightness dimming effects, and 
are also somewhat less dependent on extinction since ob- 
servations mainly probe rest-frame optical instead of UV. 
It is therefore of high interest to extend SN searches to 
larger areas, or number of observed epochs, to increase 
the sample and hence minimize the statistical part of the 
errors. A caveat is, however, that CC SNe are relatively 
faint, and it will be challenging to derive a sufficiently 
large, as well as complete, sample at redshifts exceeding 
z ~ 1 — 2, until the next generation of telescopes are 
operational. 

The results from type la rates clearly show an increase 
in the rate toz~ 1, thereafter the rate flattens and shows 
a suggested decrease in the redshift range 1 ;$ z ;$ 1.6. De- 
riving a SNR with this shape from a SFR scenario that is 
consistent with recent observations of the UV-luminosity 
density, requires that there is a significant delay between 
the star formation and the subsequent explosion of the la 
SNe. In Strolger et al. (2004), we showed that a 'narrow' 
Gaussian distribution of delay times best fits our obser- 
vations. However, other functional forms for the delay 
time may also be consistent with observations as long as 
the characteristic time delay, r, is £ 2 Gyr. Here we 
have shown that the 'narrow' Gaussian distribution is best- 
fitting also when we enlarge our sample of observed rates 
by rates taken from the literature. 

To further constrain the characteristic time delay t, as 
well as the functional forms of the delay time, it is neces- 
sary to determine the shape of the SFR to greater detail 
than known today. Especially, an increased redshift cov- 
erage will make it possible to set firmer constraints on r. 
This is shown in Figure 2, where models diverge signifi- 
cantly at higher redshift z £> 2. Alternatively, if studies 
of SN la can reveal the mechanisms that lead to the ex- 
plosion, i.e., determining the distribution of delay times, 
then it will be possible to set firmer constraints on the 
SFR using the observed SNR. 

7. CONCLUSIONS AND SUMMARY 

Using a unique sample of high redshift supernovae ob- 
tained with the ACS on-board HST, we have calculated 
the rates of core collapse SNe toz~ 0.7 and type la SNe 
to z ~ 1.6. These redshifts are significantly higher than 
those where rates have previously been measured. We find 
a core collapse SN rate that is a factor ~ 1.6 higher at 
< z >^ 0.7 than at < z >~ 0.3. Compared to local rates, 



we find an increase in the rate by a factor ~ 7 over the 
redshift range < z < 0.7. Relating the SN rate with the 
star formation rate, we find that this increase is consis- 
tent with the evolution of the cosmic star formation rate 
measured from extinction corrected galaxy UV-luminosity 
densities. The rates we derive are also consistent with the 
star formation rate estimated from mid- and far-IR lumi- 
nosity densities, even though the data do not fit as well as 
with the UV-luminosity density. A better fit with the IR 
derived star formation rate would be achieved in a scenario 
where the dust extinction increases with redshift. 

These first measurements of the SFR based on CC SNe 
are affected by both statistical uncertainties and possible 
systematic errors, as extensively discussed above. How- 
ever, most of these errors are independent of redshift. 
Therefore, the relative increase in the SFR with redshift 
is robust. 

The type la SN rate shows a significant increase to z ~ 1, 
compared to measurements at z < 0.5 taken from the lit- 
erature. In the redshift range 1 ;$ z £ 1.6, we find a SNR 
that is consistent with decreasing. When fitting observed 
rates to those derived from observed star formation rates 
combined with different assumptions on the SN la delay 
time distributions, we find that the observed shape of the 
la rate is most consistent with models assuming a signifi- 
cant delay time between the epoch of star formation and 
the explosion of the la SN. We find characteristic delay 
times t £ 2 Gyr are consistent with data. This provides 
further support to arguments favoring single-degenerate 
scenarios for the progenitors of SNe la (Livio 2001). 

Fitting observations to models also provide estimates 
of the fraction of formed stars that eventually explode as 
SN la. Assuming a Salpeter IMF and a progenitor main 
sequence mass range 3M Q < M < SMq, we find that 
5 — 7% of the white dwarfs formed from these progeni- 
tor stars subsequently explode as SNe la. Furthermore, 
if we assume that all stars with main sequence masses 
O.8M < M < 8M Q become white dwarfs, then 0.6-0.8% 
of these will turn into SNe. Investigating the intrinsic ra- 
tio of exploding CC to la SNe results in a ratio ~ 3, which 
is fairly constant up to z ~ 0.7. At higher redshifts the 
ratio increases, reflecting the fact that a long delay time 
for la progenitors shifts the rate towards lower redshift 
compared to the SFR and therefore also the CC SNR. 

To set firm constraints on the SFR using core collapse 
SNe, as well as determining the distribution of delay times 
for SNe la progenitors, it is of great importance not only 
to increase the sample size (i.e. decrease statistical errors), 
but also to better understand SN characteristics. Never- 
theless, using available information on SN statistics and 
characteristics, combined with the unique sample of high 
redshift SNe provided by the GOODS, we have already 
shown that using high redshift supernova rates can pro- 
vide important information about the properties and the 
evolution of stars and star formation in the Universe at 
earlier cosmic epochs. 
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Fig. 1. — Core collapse supernova rates from GOODS are shown as filled circles. Vertical error bars show statistical errors, while horizontal 
error bars show bin size. Star formation model Ml (solid black line) is derived from the data in Giavalisco et al. (2004b), and uses similar 
assumptions about dust extinction as applied here. The star formation rate derived from analysis of various mid- and far- infrared datascts 
by Chary & Elbaz (2001) is shown with the solid gray line. Also shown are star formation model M2 (dashed line) which is also derived 
from the Giavalisco et al. data and a star formation rate model given in Lanzetta et al. (2002) (dotted line). The latter two models are not 
corrected for dust extinction. For a comparison between these and our results, we show as open circles the rates derived assuming no dust 
extinction. The local rate from Cappellaro et al. (1999) is shown as a filled square. 
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Fig. 2. — Type la SNRs from GOODS at < z >=0.40, 0.80, 1.20, and 1.60 arc shown as filled circles. Also plotted are rates at z ~ 0.01 from 
Cappellaro et al. (1999)(fillcd square), at z ~ 0.1 from Hardin ct al. (2000)(open star), at z ~ 0.11 from Strolgcr ct al. (2003a) (filled star), 
at z ~ 0.11 from Rciss (2000)(opcn triangle), at z ~ 0.38 from Pain et al. (1996)(opcn square), at 2 ~ 0.46 from Tonry ct al. (2003)(opcn 
circle), and at z ~ 0.55 from Pain et al. (2002)(fillcd triangle). Vertical error-bars on the GOODS rates represent statistical errors, while 
horizontal error-bars represent bin size. See text for discussion on systematic errors. The figure also shows predicted la rates based on three 
different models for the delay time distribution of SN la progenitors. 
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Fig. 3. — Intrinsic ratio of exploding core collapse to la SNe as a function of redshift. Ratio is shown for SFR model Ml using three different 
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Table 1 
Supernova characteristics. 



Typo 


Mb 


a 


/ 




Type la 






Ia-'bright' 


-19.6 


0.30 


0.20 


Ia-'normal' 


-19.3 


0.45 


0.64 


Ia-'faint' 


-17.8 


0.50 


0.16 




Core collap 


se 




Ibc 


-17.28 


0.74 


0.225 


Ibc-'bright' 


-19.93 


0.33 


0.005 


IIP 


-16.66 


1.12 


0.300 


IIP-'faint' 


-14.39 


1.00 


0.150 


IIL 


-17.22 


0.38 


0.295 


IIL-'bright 


-18.94 


0.51 


0.005 


Iln 


-18.82 


0.92 


0.020 



Note. — Columns are: (1) Supernova 
type; (2) Absolute peak B band mag- 
nitude assuming Hq—70 km s~ Mpc - 
and a cosmology with Q,m — 0.3 and 
fl\ — 0.7; (3) Dispersion in peak magni- 
tude; (4) Intrinsic fractions of subtypes 
within the two groups. 
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Table 2 
Supernova rates 



Rcdshift: 0.2 < z < 0.6 0.6 < z < 1.0 1.0 < z < 1.4 1.4 < z < 1.8 

Twe la Q gg +0.34 +1.54 1 „ +0.44 +0.75 1 _ +0.47 +0.32 „ 44 +0.32 +0.14 

lype la u.vv _ ,27 -0.25 1 -° l -0.25 -0.53 A,xu -0.26 -0.44 u.<±<± _o.25 -0.11 

Number of SNe 3 14 6 2 

Redshift: 0.1 < z < 0.5 0.5 < z < 0.9 

Core collapse 2.51 to.jl ti'.se 3 - 96 il'.oe -2'.60 
Rates without extinction correction 

. ._ +0.60 +0.42 , „„ +0.53 +0.65 

1 - ' -0.46 -1.12 1 ' DD -0.36 -0.93 

Number of SNe 6 10 



Note. — Rates are given in units yr~ Mpc~ 10~ ^705 assuming a cosmology with Qm — 
0.3 and Cl\ — 0.7. First quoted errors are statistical and represent 68.3% confidence intervals. 
Second errors are systematic and include the possibility that all SNe with uncertainty in type 
determination have been misclassificd, as well as other possible sources discussed in §5.1.1 and 
§5.2.1. Note that the latter errors are non-Gaussian. 



